SCALING LAWS FOR NON-EUCLIDEAN PLATES AND THE W 2 > 2 
ISOMETRIC IMMERSIONS OF RIEMANNIAN METRICS 

MARTA LEWICKA AND REZA PAKZAD 

Abstract. This paper concerns the elastic structures which exhibit non-zero strain at free equi- 
libria. Many growing tissues (leaves, flowers or marine invertebrates) attain complicated config- 
urations during their free growth. Our study departs from the 3d incompatible elasticity theory, 
conjectured to explain the mechanism for the spontaneous formation of non-Euclidean metrics. 

Recall that a smooth Riemannian metric on a simply connected domain can be realized as the 
pull-back metric of an orientation preserving deformation if and only if the associated Riemann 
curvature tensor vanishes identically. When this condition fails, one seeks a deformation yielding 
the closest metric realization. We set up a variational formulation of this problem by introducing 
the non-Euclidean version of the nonlinear elasticity functional, and establish its F-convergence 
under the proper scaling. As a corollary, we obtain new necessary and sufficient conditions for 
existence of a W 2 ' 2 isometric immersion of a given 2d metric into K 3 . 
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1. Introduction 

Recently, there has been a growing interest in the study of flat thin sheets which assume 
non-trivial configuration in the absence of exterior forces or impose boundary conditions. This 
phenomenon has been observed in different contexts: growing leaves, torn plastic sheets and 
specifically engineered polymer gels [14 . The study of wavy patterns along the edges of a torn 
plastic sheet or the ruffled edges of leaves suggest that the sheet endeavors to reach a non-attainable 
equilibrium and hence necessarily assumes a non-zero stress rest configuration. 

In this paper, we attempt to give a possible mathematical foundation of these phenomena, 
in the context of the nonlinear theory of elasticity. The basic model, called "three dimensional 
incompatible elasticity" |4;, follows the findings of an experiment described in p3]. The experi- 
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Figure 1.1. The experimental system and the obtained structures of sheets with 
radially symmetric target metrics. Reprinted from [J3] with permission from 
A A AS. 



ment (see Figure 1.1) consists in fabricating programmed flat disks of gels having a non-constant 
monomer concentration which induces a "differential shrinking factor" . The disk is then activated 
in a temperature raised above a critical threshold, whereas the gel shrinks with a factor propor- 
tional to its concentration and the distances between the points on the surface are changed. This 
defines a new target metric on the disk, inducing hence a 3d configuration in the initially planar 
plate. One of the most remarkable features of this deformation is the onset of some "transversal" 
oscillations (wavy patterns). 

Trying to understand the above phenomena in the context of nonlinear elasticity theory, it has 
been postulated |14[ [1] that the 3d elastic body seeks to realize a configuration with a prescribed 
pull-back metric g. In this line, a 3d energy functional was introduced in jj], measuring the L 2 
distance of the realized pull-back metric of the given deformation from g. 

Unfortunately, this functional is not suitable from a variational point of view, due to the exis- 
tence of Lipschitz deformations with zero energy level [7j . These deformations are however neither 
orientation preserving nor reversing in any neighborhood of a point where the Riemann curvature 
of the metric g does not vanish (i.e. when the metric is non-Euclidean). In order to overcome this 
shortcoming, here we introduce a modified energy I{u) which measures, in an average pointwise 
manner, how far a given deformation u is from being an orientation preserving realization of the 
prescribed metric. An immediate consequence is that for non-Euclidean g, the infimum of / (in 
absence of any forces or boundary conditions) is strictly positive, which points to the existence of 
non-zero strain at free equilibria. 

Several interesting questions arise in the study of the proposed energy functional. A first one 
is to determine the scaling of the infimum energy in terms of the vanishing thickness of a sheet. 
Another is to determine limiting zero-thickness theories under obtained scaling laws. The natural 
analytical tool in this regard is that of T convergence, in the context of Calculus of Variations. 

In this paper, we consider a first case where the prescribed metric is given by a tangential 
Riemannian metric \g a p] on the 2d mid-plate, and is independent of the thickness variable. The 
3d metric g is set-up such that no stretching may happen in the direction normal to the sheet in 
order to realize the metric. Consequently, if \g a p\ has non-zero Gaussian curvature, then such g 
is non-Euclidean. We further observe a correspondence between the scaling law for the infimum 
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energy of the thin sheet in terms of the thickness, and the immersability of [g a p] into (Theorem 



2.6). This result relates to a longstanding problem in differential geometry, depending heavily on 
the regularity of the immersion |15| [9"]. In our context, we deal with W 2 ' 2 immersions not studied 
previously. We also derive the T-limit of the rescaled energies, expressed by a curvature functional 



on the space of all W 2,2 realizations of [g a p] m ^ 3 (Theorems 2.4 and 2.5). 

To put our results in another context, recall the seminal work of Friesecke, James and Miiller 
[5], where the nonlinear bending theory of plates (due to Kirchhoff) was derived as the r-limit 
of the classic theory of nonlinear elasticity, under the assumption that the later energy per unit 
thickness h scales like h 2 . From a mathematical point of view, the present paper provides the 
non-Euclidean version of the same results under the same scaling law, and the 2d limit theory we 
obtain is the natural non-Euclidean generalization of the Kirchhoff model. Contrary to the classic 
case, in our context the scaling law is the unique natural scaling of the energy for the free thin 
sheet with the associated prescribed metric. 

As a major ingredient of proofs, we also give a generalization of the geometric rigidity estimate 
[5] to the non-Euclidean setting (Theorem 2.3). We estimate the average I? oscillations of the 



deformation gradient from a fixed matrix in terms of the 3d non-Euclidean energy / and certain 
geometric parameters of the 3d domain. The main difference is an extra term of the bound, 
depending on the derivatives of the prescribed metric g and hence vanishing when g is Euclidean 
as in [5] . 

Acknowledgments. We are grateful to Stefan Miiller for a significant shortening of the original 
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2. Overview of the main results 

Consider an open, bounded, Lipschitz domain IA C M n , with a given smooth Riemannian metric 
g = [gij]. The matrix field g : U — * M nxn is therefore symmetric and strictly positive definite up 
to the boundary dli. Let A = ^J~g be the unique symmetric positive definite square root of g and 
define, for all x G U: 

(2.1) T(x) = \rA{x); R G SO(n)\, 

where SO{n) stands for the special orthogonal group of rotations in W 1 . By polar decomposition 
theorem, it easily follows that u is an orientation preserving realization of g: 

(Vu) T Vu = g and det Vu > a.e. in U 

if and only if: 

Vu(x) G ^(x) a.e. in U. 
Motivated by this observation, we define: 

(2.2) I(u) = [ dist 2 (Vu(x),F(x)) dx \/u G W 1,2 (U,M n ). 

Ju 



Notice that when g = Id then the above functional becomes I(u) = J dist 2 (Vu, SO(n)) which is 
a standard quadratic nonlinear elasticity energy, obeying the frame invariance. 
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Remark 2.1. For a deformation u : hi — ► W 1 one could define the energy as the difference 
between its pull-back metric on hi and g: 



IstAu) = [ \{Vu) T Vu-g\ 2 dx. 
JU 



However, such "stretching" functional is not appropriate from the variational point of view, for the 
following reason. It is known that there always exists u G W 1 ' co (U,W l ) such that I s t r (u) = 0. On 
the other hand [7] , if the Riemann curvature tensor R associated to g does not vanish identically, 
say Rijki(x) ^ for some x £U, then u must have a 'folding structure' around x; the realization 
u cannot be orientation preserving (or reversing) in any open neighborhood of x. 



In view of the above remark, our first observation concerns the energy (2.2 ) in case of R 0. 



Theorem 2.2. If the Riemann curvature tensor Rijki ^ 0, then inf |/(u); u G W 1,2 (JA, M n )| > 0. 

In case g = Id, the infimum as above is naturally and is attained only by the rigid motions. 
In |5j, the authors proved an optimal estimate of the deviation in W 1,2 of a deformation u from 
rigid motions in terms of I(u). In section [4] we give a generalization of such quantitative rigidity 
estimate to our setting. Since there is no realization of I(u) = if the Riemann curvature is non- 
zero, we choose to estimate the deviation of the deformation from a linear map at the expense of 
an extra term, proportional to the gradient of g. 

Theorem 2.3. For every u G W 1,2 (U,M. n ) there exists Q £ R nxn such that: 

\Vu{x)-Q\ 2 dx<c(j dist 2 (Vn,J^(x)) dx + || Vg\\ 2 L oo (diam U) 2 \U\^j , 

where the constant C depends on ||<7||l°°> a nd on the domainlA. The dependence on 

IA is uniform for a family of domains which are bilipschitz equivalent with controlled Lipschitz 
constants. 

We shall consider a class of more general 3d non-Euclidean elasticity functionals: 

I w {u) = / W(x,Vu(x)) dx, 
Ju 

where the inhomogeneous stored energy density W : li x M nxri — > M + satisfies the following 
assumptions of frame invariance, normalization, growth and regularity: 

(i) W{x,RF) = W(x,F) for all R G SO(n), 

(ii) W(x,A(x)) = 0, 

(iii) W(x,F) > c dist 2 {F,T(x)), with some uniform constant c > 0, 

(iv) W has regularity C 2 in some neighborhood of the set {(x, F); x £U,F G ^(x)}. 

The properties (i) - (iii) are assumed to hold for all x G IA and all F G M. nxn . In case the Riemann 



curvature tensor of g does not vanish, by Theorem 2.2 the infimum of I\y is positive, in which 
case Iw is called a three dimensional incompatible elasticity functional. 

We consider thin 3d plates of the form: 

n h = n x {-h/2, h/2) C M 3 , 0<h«l, 

with a given mid-plate Q an open bounded subset of M 2 . In accordance with [T3] , we assume that 
the metric g on O' 1 has the form: 





(2.3) 9(x',x 3 ) 
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where [g a (3\ is a smooth metric on ft, defined up to the boundary. In particular, g does not depend 
on the thin variable x%. Accordingly, we shall assume that the energy density W does not depend 
on X3: 

(v) W(x, F) = W(x', F), for all x G U and all F G R nxn . 
Define now the rescaled energy functionals: 



I h ( u ) = l [ W(x,Vu(x)) dx VueW 1 ' 2 ^ 



where the energy well F(x) = ^(x') = SO(3)A(x) is given through the unique positive definite 
square root A = ^Jg~ of the form: 



A(x',x 3 ) 



A a(3 {x') 
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Vx' 6 fi, X3 G (—h/2, h/2). 



By an easy direct calculation, one notices that the Riemann curvature tensor Rijki = 0, of g in 
£l h if and only if the Gaussian curvature of the 2d metric K[ 9q(3 ] = 0. Hence, by Theorem 



2.2 



inf I h > for all h if this condition is violated. A natural question is now to investigate the 
behavior of the sequence inf/' 1 as h — > 0. We first obtain (in section |5| the following lower bound 
and compactness result: 

Theorem 2.4. Assume that a given sequence of deformations u h G W 1,2 (£~l h , M. 3 ) satisfies: 
(2.4) I h (u h ) < Ch 2 , 

where C > is a uniform constant. Then, for some sequence of constants c h <E M 3 , the following 
holds for the renormalized deformations y h (x' ,x 3 ) = u h (x' ,hx 3 ) — c h G P^ 1,2 (f2 1 ,R 3 ): 

(i) y h converge, up to a subsequence, in W 1,2 (f2 1 , M 3 ) toy(x',x 3 ) = y(x') andy G W 2,2 (fl, IR 3 ). 

(ii) The matrix field Q(x') with columns Q(x') 
x' G Here: 

d\y x d 2 y 



diy(x'),d 2 y(x'),n(x') G ^(x 1 ), for a.e. 



(2.5) 



n 



\diy x d 2 y\ 



is the (well defined) normal to the image surface y(£l). Consequently, y realizes the mid- 
plate metric: (Vy) T Vy = [g a p\. 
(iii) Define the following quadratic forms: 



Q 3 (x')(F) = X7 2 W(x', -)\ A(xl) {F, F), Q 2 {x'){F tan ) = mm{Q 3 (x')(F); F tan 
Then we have the lower bound: 



Ftan } • 



Ynn\n{^I h (u h ) > 
h^o h 2 v ' ~ 24 



l -j^Q 2 { x ')[A- a l{Vy) T Vn) dx' . 



We further prove that the lower bound in (iii) above is optimal, in the following sense. Let 
y G W 2 ' 2 (£l, M 3 ) be a Sobolev regular isometric immersion of the given mid -plate metric, that is 
(X7y) T X7y = [g a f3\- The normal vector n G W l,2 (Q, M 3 ) is then given by (2.5) and it is well defined 
because \d±y x d 2 y\ = (det g) 1 ^ 2 > 0. 
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Theorem 2.5. For every isometric immersion y E W > ) of g, there exists a seq uence 



of "recovery" deformations u h E W 1,2 (Q h , M 3 ) such that the assertion (i) of Theorem 2.4 
together with: 

(2.6) lim ^l\u h ) = 1 ! Q 2 {x>) (A- a l{Vy) T Vn) Ax'. 



hold, 



A corollary of Theorems 2.4 and 2.5 proved in section [8] provides a necessary and sufficient 
condition for the existence of W 2,2 isometric immersions of (O, [<? Q/ 3]): 

Theorem 2.6. Let [g a /3] be a smooth metric on the midplate $7 C M 2 . Then: 

(i) [Qafi] has an isometric immersion y £ VF 2,2 (f2, M 3 ) if and only if ^ ml I h < C, for a 
uniform constant C . 

(ii) [g a p\ has an isometric immersion y £ W 2 ' 2 (Q, M 2 ) (or, equivalently, the Gaussian curva- 
ture ^[g al3 \ = 0) if and only i/lim^— >o m f I h = 0. 

(hi) If the Gaussian curvature K\g a/3 ] ^ in t/ien 7^ inf I h > c > 0. 

The existence (or lack of thereof) of local or global isometric immersions of a given 2d Riemann- 
ian manifold into M 3 is a longstanding problem in differential geometry, its main feature being 
finding the optimal regularity. By a classical result of Kuiper |15j . a C 1 isometric embedding into 
M 3 can be obtained by means of convex integration (see also [7])- This regularity is far from W 2,2 , 
where information about the second derivatives is also available. On the other hand, a smooth 
isometry exists for some special cases, e.g. for smooth metrics with uniformly positive or negative 
Gaussian curvatures on bounded domains in M 2 (see [9], Theorems 9.0.1 and 10.0.2). Counterex- 
amples to such theories are largely unexplored. By [12] , there exists an analytic metric [g a p\ with 
nonnegative Gaussian curvature on 2d sphere, with no C 3 isometric embedding. However such 
metric always admits a C ' embedding (see [8] and [TO] ) . For a related example see also [22]. 



Finally, notice that Theorems 2.4 and |2.5| can be summarized using the language of T-convergence 
PJ- Recall that a sequence of functionals T : X — > R defined on a metric space X, is said to 
T-converge, as h — > 0, to T : X — > M provided that the following two conditions hold: 

(i) For any converging sequence {x h } in X: 

T\ lima;' 1 ) < liming (x fc ). 
\h-+o ) h^O 

(ii) For every x £ X, there exists a sequence {x h } converging to x and such that: 

T(x) = \imF h (x h ). 

Corollary 2.7. The sequence of functionals T : W 1)2 {Q}, M 3 ) — > M, given by: 

f h (y(x)) = ±I h (y(x',hx 3 )) 

Y -converges, as h — > 0, to: 

I +00 otherwise. 
Consequently, the (global) approximate minimizers of T h converge to a global minimizer of T . 
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3. The non- Euclidean elasticity functional - a proof of Theorem 12.21 

In the sequel, we will need some differential geometry notation. We shall write \g\ = det g and 
g^ 1 = [g tJ ']. The Christoffel symbols are then given, using the Einstein summation, as: 

T^ = ^g km (d i g kj + d j9ik -d k g ij ). 

By V we denote the covariant gradient of a scalar/vector field or a differential form, while by V g 
we denote the contravariant gradient. The covariant divergence of a vector field u can be written 
as: 

div g u = (Viu) 1 = —=di(^/\g\u l ) 
V\9\ 

and the scalar product of two vector fields (that is of two (1,0) contravariant tensors) has the 
form: (u, v) g = tfgijV 3 . We shall often use the Laplace-Beltrami operator A g of scalar fields /: 

A g f = div 5 (V 9 /) = *9i{y/\g\<f3djf). 

v m 



By R = [Rijki] we mean the ((0,4) covariant) Riemann curvature tensor, and by Ric g the ((0,2) 
covariant) Ricci curvature tensor. 

Towards the proof of Theorem |2.2[ we first derive an auxiliary result, which is somewhat 
standard in differential geometry (see e.g. |23[ [2]). 

Lemma 3.1. Let u E W 1 ' (14, W 1 ) satisfy Vu(x) E J-(x) for a. a. x E U. Then u is smooth and 
R = 0. 

Proof. Write u = (it 1 , . . . , u n ) and notice that in view of the assumption, each u l E W 1,co . More- 
over: 

det Vu = \f\g~\i cof Vu = y^H(Vu)5 _1 . 
Recall that for a matrix F E M nxn , cof F denotes the matrix of cofactors of F, that is (cof F)ij = 
(— 1)* +J det Fjj, where Fij E ]R( n_1 ) x ( n ~ 1 ) is obtained from F by deleting its ith row and jth 
column. Since div(cof Vu) = (the divergence of the cofactor matrix is always taken row- wise), 
the Laplace-Beltrami operator of each component u m is zero: 

AgU m = 0, 

and therefore we conclude that u m E C°°. The second statement follows immediately since u : 
IA — > M. n is a smooth isometric embedding of (U, g) into the Euclidean space W 1 . ■ 

Remark 3.2. For the convenience of the reader we now give a simple argument proving that the 
existence of a smooth u as in Lemma 3.1 implies that the Ricci curvature tensor Ric g = 0. Recall 
that when n = 3 (which is the dimension relevant to our main results), Ric g = if and only if 

r = o$s\. 

We shall first deduce that the second g-covariant derivative of each scalar field u m vanishes. 
Since the vectors {d s u}™ =1 form a basis of W 1 (Vu being invertible), it is enough to consider the 
following linear combination of components of V 2 u m (for fixed covariant indices 

n 

V^Vu™), • d s u m = (df jU - T%d k u)d s u 

m=l 

1 

2 

dfjud s u - -g kl (djg is + dig js - d s gij) = 0. 



dfjud s u - -g kl (djgu + digji - digij)d k ud s u 
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Hence V 2 u m = and since |V 2 u m |^ = IV 2 ,?/™) 2 ,, we also see that the second contravariant gradient 
of each component u m vanishes: 



(3.1) V z g u m = 0. 

On the other hand, (Vu) T (Vu)(7 _1 = Id, so (Vu)5~ 1 (Vm) t = Id which means precisely that 
diU m g^djU m = 1. Therefore: 

(3.2) |V 5 -u m |2 = \Vu m \ 2 g = 1. 



Applying now (3.1) and (3.2) in the following Bochner's formula (see e.g. |21j): 

\A g \V 9 f\ 2 = (V g A g f,V a f) 9 + |V 2 /| 2 + Ric g (V g f,V g f), 

where we take the scalar field / = u m , we obtain: 

Ric g (V g u m ,V g u m ) = 0. 

Since {V 5 tt m }^ =1 form a basis of M. n and Ric g is a symmetric bilinear form, we conclude that 
Ric g = 0. ■ 

We now prove two further auxiliary results. 

Lemma 3.3. There is a constant M > 0, depending only on \\g\\L°° and such that for every 
u £ W ,2 {U, K n ) there exists a truncation u G W l,2 iU,W l ) with the properties: 

||Vu||i,« < M, \\Vu - Vu||| 2(w s < 4J(u) and I(u) < WI(u). 

Proof. Use the approximation result of Proposition A.l. in [5 to obtain the truncation u = u x , 
for A > having the property that if a matrix F G M nxn satisfies \F\ > A then: 

\F\ <2dist 2 (F,F{x)) VxeU. 

Then ||Vu A || L oo < CA := M and further: 

\\Vu-Vu x \\ 2 L 2 < \Vu\ 2 <4 dist 2 (Vu,F(x)) dx <4I(u). 

i{|V«|>A} i{|V«|>A} 
The last inequality of the lemma follows from the above by triangle inequality. ■ 

Lemma 3.4. Let u G W 1,oa {U, W l ) and define vector field w whose each component w m satisfies: 

A g w m = inU, w m = u m ondU. 

Then ||V(n — t^OH 2 ^^ — d( u )> where the constant C depends only on the coercivity constant of 
g and (in a nondecreasing manner) on ||Vii||i,°°. 



Proof. The unique solvability of the elliptic problem in the statement follows by the usual Lax- 
Milgram and compactness arguments. Further, the correction z = u — w G W ' 2 (jU, W 1 ) satisfies: 

/ 9 ii ^\g\d i z m d j 4>= [ g ij V\9\d iU m d j( l>- [ V0(cofVn) m _ th 
Ju Ju Ju 



row 
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1 2 

for all <j) G W ' (Li '). Indeed, the last term above equals to 0, since the row- wise divergence of the 
cofactor matrix of Vu is 0, in view of u being Lipschitz continuous. Use now <f> = z m to obtain: 



m\2 



u 



m\2 



U 



(3.3) 



/ 
Ju 



in 

< \\Vz m \ 



d jZ m (g ij ^/\g\diu m - (cof Vu) mj 



J 

Ju 



\9\(Vu)g 



-i 



cof Vu 



1/2 



<C M \\Vz m \\ L HuAu) ll \ 



which proves the lemma. 

In order to deduce the last bound in (3.3), consider the function f(F) = \Z\g\Fg~ 1 — cof F, 
which is locally Lipschitz continuous, uniformly in x G Vi. Clearly, when F G J~(x) then F = RA 
for some R G SO(n), and so: cof F = cof (RA) = (det A) RA' 1 = y /\g\(RA)g- 1 , implying: 
f(F) = 0. Hence: 

|/(Vn(x))| 2 < C 2 M dist 2 (Vu(x),F(x)), 

where Cm stands for the Lipschitz constant of / on a sufficiently large ball, whose radius is 
determined by the bound M = ||Vu||l<». ■ 



2.2\ We argue by contradiction, assuming that for some sequence of 



Proof of Theorem 

deformations u n G W 1,2 (U, W 1 ), there holds limn^oo I(u n ) = 0. By Lemma 3.3, replacing u n by 
u n , we may also and without loss of generality have u n G W ' 0O (U,M. n ) and ||Vu n ||Loo < M. 

Clearly, the uniform boundedness of Vu n implies, via the Poincare inequality, after a modifica- 
tion by a constant and after passing to a subsequence if necessary: 

(3.4) limii n = u weakly in W 1,2 (U). 



Consider the splitting u n = w n + z n as in Lemma 3.4 By the Poincare inequality, Lemma 3.4 

limz„ = strongly in W 1,2 (U). 



implies that the sequence z n G W ' 2 (U) converges to 0: 



In view of the convergence in (3.4), the sequence w n must be uniformly bounded in W > 2 (U), 
and hence by the local elliptic estimates for the Laplace-Beltrami operator, each A g -harmonic 
component w™ is locally uniformly bounded in a higher Sobolev norm: 

W' CCU 3Cw \\ w n\\w 2 > 2 (W) < C'w'll u Cllw 1 . 2 (w) < C. 

Consequently, w n converge to u strongly in W^ 2 (U) and recalling that I(u n ) converge to 0, we 
finally obtain: 

I(u) = 0. 

Therefore Vu G J~(x) for a. a. x GU, which achieves the desired contradiction with the assumption 
R ^ 0, by Lemma [O ■ 
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4. A GEOMETRIC RIGIDITY ESTIMATE FOR RlEMANNIAN METRICS - A PROOF OF THEOREM 12.31 

Recall that according to the basic rigidity estimate [5], for every v € W Al,2 (V,M n ) defined on an 
open, bounded set V C W 1 , there exists R £ SO(n) such that 



n)). 



(4.1) I \Vv-R\ 2 <C v [ dist 2 {Vv, SO{ 

Jv Jv 

The constant Cy depends only on the domain V and it is uniform for a family of domains which 
are bilipschitz equivalent with controlled Lipschitz constants. 



A proof of Theorem 2.3, For some xo £ U denote Aq = A(xq) and apply (4.1) to the vector 
field v(y) = u(Aq 1 ij) £ W L ' 2 (AoU, M n ). After change of variables we obtain: 

3R G SO{n) [ \(Vu)Aq 1 - R\ 2 < C Aq u I dist 2 ((Vu),4 \ SO(n)). 
Ju Ju 

Since the set AqU is a bilipschitz image of U, the constant Ca u has a uniform bound C depending 

on H^oll, II^q 1 !! an d ZY- Further: 



f \Vu - RA \ 2 < C\\A \\ 4 [ dist 2 (Vu, SO{n)A ) 
Ju Ju 

<C\\A \\ 4 ( [ dist 2 {Vu,F{x)) + f \A{x)-A 
\Ju Ju 

<C||ffll!°° (j dist 2 (Vn,J^(x)) dx + C||Vc/|||oo(diam^) 2 |W| 

which proves the claim. ■ 

We now derive a crucial approximation result, as in Theorem 10 [6] (see also Lemma 8.1 [18 J. 
Lemma 4.1. There exist matrix fields Q h G W l,2 (VL, M 3x3 ) such that: 

(4.2) \[ \Vu h (x) -Q h (x')\ 2 dx <C(h 2 + I h (u h )), 

(4.3) f \VQ h \ 2 <C(l + h- 2 I h (u h )), 

Jn 

with constant C independent of h. 

Proof. Let D x i ^ = B(x', /i)flfl be 2d curvilinear discs in ft of radius h and centered at a given x 1 . 
Let B x ' t h = D x r t h x (—h/2,h/2) be the corresponding 3d cylinders. On each B x i^ use Theorem 
12.31 to obtain: 



f \Vu h -Q x , :h \ 2 <CU dist 2 (V u h ,F(z)) dz + h 2 \B x , !h \ J 
(44) jB x',h V B x',h ) 

<C f h 2 + dist 2 (ViA F{z)) dz, 

with a universal constant C in the right hand side above, depending only on the metric g and the 
Lipschitz constant of d£l. 

Consider now the family of mollifiers r\ x i : Q — ► R, parametrized by i' E fl and given by: 

d(\z'-x'\/h) 



r)x'(z') 



hf n 9(\y'-x'\/h) dy'' 
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where 9 E C£°([0, 1)) is a nonnegative cut-off function, equal to a nonzero constant in a neighbor- 
hood of 0. Then rj x '(z') = for all z 1 D x ^ and: 



/ 



Vx> = h X , H^x'IIl 00 < C/i 11 , 1 1 Vx'Tjx' 1 1 < Ch 



-3 



Define the approximation Q n G VF 1,2 (0, M 3x3 ): 



Q h (x') = / ij^z'jViiV) dz. 



n h 



By (4.4) , we obtain the following pointwise estimates, for every x' £ £1 

2 



\Q h (x) - Q x >,h\ 2 < U Vx'(z') (Vu h (z) - Q x < t hj dz 

< [ \r]x>(z')\ 2 dz ■ [ \Vu h - Q x > h \ 2 < Ch~ 3 [ h 2 + dist 2 (V u h ,F(z)) dz, 

\VQ\x')\ 2 = U(V xl r, x ,(z'))Vu h (z) dz)' 

= ([ (V x >Vx>(z'))(Vu h (z)-Q x , th ) dz) <[ \V x ,r, x ,(z')\ 2 dz ■ [ \Vu h -Q x ,, h \ 2 
\Jn h v ' J Jn h Jn h 

< Ch~ 5 [ h 2 + dist 2 (W\ F(z)) dz. 
Applying the same estimate on doubled balls B x i^h we arrive at: 

/ \Vu h (x)-Q h (x')\ 2 dx<C[ [ \Vu h (z)-Q x ,, h \ 2 dz+ [ \Q x ,, h -Q 

<C I h 2 + dist 2 (Vu h , F{z)) dz, 

J B x , ah 

j \VQ h \ 2 < Ch' 3 j h 2 + dist 2 (Vu h , F(z)) dz. 

Consider a finite covering O = (J D x /^ whose intersection number is independent of h (as it 
depends only on the Lipschitz constant of d£l). Summing the abov e bo unds and applying the 
uniform lower bound W(x, F) > c dist 2 (-F, J~(x)) directly yields (4.2) and (4.3). ■ 



h (z')\ 2 dz 



5. Compactness and the lower bound on rescaled energies - a proof of Theorem 

El 



1. From (4.2), (4.3) in Lemma 4.1 and the a ssumption on the energy scaling, it follows directly 
that the sequence Q h , obtained in Lemma 4.1 is bounded in W 1,2 (0,, IR 3x3 ). Hence, Q h converges 



weakly in this space, to some matrix field Q and: 

/ \Vu h (x',hx 3 ) -Q(x')\ 2 dx< [ \Q h -Q\ 2 + h- 1 I \Vu h (x) - Q h (x')\ dx, 
Jn 1 Jo, Jn h 



converges to by (4.2). Therefore we obtain the following convergence of the matrix field with 
given columns: 

lim \d 1 y h (x),d 2 y h (x),h- 1 d 3 y h (x)\ = Q(x) in L 2 (U\ R 3x3 ). 

h— >0 L 
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We immediately conclude that ||53?/ /l ||/ / 2(Qi) converges to 0. 

Now, setting c h = u h (x', hx 3 ) dx, by means of the Poincare inequality there follows the 
assertion (i) of the Theorem. The higher regularity of y can be deduced from Vy G W 1,2 , in view 
of the established W 1,2 regularity of the limiting approximant Q. 

To prove (ii), notice that by (4.2), (4.4) and the lower bound on W: 

(5.1) f dist 2 (Q h ,T(x')) < Ch 2 . 



Hence Q(x') G Fix) a.e. in Q and consequently d a y ■ dpy = g a p. To see that the last column of 
the matrix field Q coincides with the unit normal to the image surface: Qe 3 = n, write Q = RA 
for some R G SO{2>) and notice that: 

diy x d 2 y = (RAei) x (RAe 2 ) = R{{Aei) x (Ae 2 )) = cRAe 3 = cQe 3 , 

where c = |<9iy x d 2 y\ = \ (Aei) x (Ae2)| = detA > 0, by the form of the matrix A. On the other 
hand \Qe 3 \ = \Re 3 \ = 1, so indeed there must be Qe 3 = n = [d\y x d 2 y)/\d\y x d 2 y\. 

2. We now modify the sequence Q h to retain its convergence properties and additionally get 
Q h (x') G T{x) for a.a. xefi. Define Q h G L 2 (ft,R 3 ) with: 

Kr( x ){Q h {x')) if Q h (x') G small ngbhood of JF(x) 
A(x) otherwise 



Q h (x" 



where Ttp( x ) denotes the projection onto the compact set J-{x) of its (sufficiently small) neighbor- 
hood. One can easily see that: 

[ \Q h - Q h \ 2 < C [ dist 2 (Q h (x'),F(x')) dx' < Ch 2 
Jn Jn 



by (5.1). In particular, Q h converge to Q in L 2 (f2). 

Write Q h = R h A for a matrix field R h G SO(3) and consider the rescaled strain: 

G h (x', x 3 ) = i ((R h f(x')Vu h (x', hx 3 ) - A(x')) G L 2 (U\ M 3x3 ). 

We obtain: 

/ \G h \ 2 <Ch~ 3 [ \Vu h -Q h \ 2 <Ch- 2 I h (u h ) + Ch~ 2 [ \Q h -Q h \ 2 <C. 
Jn 1 Jn h Jn 

Hence there exists a limit: 

(5.2) limG fc = G weakly in L 2 (ft\ M 3x3 ). 

h— >0 

3. Fix now a small s > and consider the difference quotients: 

f> h (x) = l \ y \ x + se 3 )- y \x)) G W 1,2 ^ 1 ,R 3 ). 
h s 

Since h~ 1 d 3 y h converges in L 2 (Q 1 ,]R 3 ) to n(x'), then also: 

1 r s 

lim f s ' h (x) = lim - f d 3 y h (x + te 3 ) dt = n(x). 

h— >0 h— >0 h Jq 

There also follows convergence of normal derivatives, strongly in L 1 (0 1 ): 

lim d 3 f s ' h (x) = lim(d 3 y h (x + se 3 ) - d 3 y h {x)) = 0, 

h— >0 h^O 



and of tangential derivatives, weakly in L (0 ): 
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1 



1 



lim d a f s ' h (x) = lim -R h (x') [G h (x + se 3 ) - GHx) )e a = -{QA~ )(x) ( G(x + se 3 ) — G(x) )e a , 

h^O h— >0 S V / S \ / 

where we have used that the L°° sequence R h = Q h A~ l converges in L 2 (Q) to QA^ 1 6 .SO^). 

Consequently, the sequence / s,?t converges, as h — » weakly in L 2 (J7 1 ) to n(x'). Equating the 
derivatives, we obtain: 



d a n(x') = l -{QA~ l ){x')(G{x + se z ) -G(x))e Q . 



Therefore: 
(5.3) 



a = 1,2. 



4. We now calculate the lower bound of the rescaled energies. To this end, define the sequence 
of characteristic functions: 

Xh = X^xen 1 ; |G h (^')l</i~ 1/2 }' 



which by (5.2) converge in L (O ) to 1. Using frame invariance and noting that in the Taylor 
expansion of the function F i— ► W(x, F) at A{x) the first two terms are 0, we obtain: 

Xh(x)W(x, X7u h (x' ,hx 3 )) dx 
Xh{x)W{x,R h {x') T Vu h {x',hx 3 )) dx 

Xh{x)W{x,A(x) + hG h (x)) dx 
1. 



1 
1 



IQ 1 

> I Xh(x) 
n 1 



V 2 W(x, .) lA{x) (G h (x), G\x)) - o(l)\G h (x)\ 



dx. 



Hence: 



liminf ^I h (u h ) > -liminf / Xh(x)Q 3 (x')(G h (x)) dx 

= \ liminf f Q 3 (x')( Xh (x)G h (x)) dx 
2 h-M) J n i V / 

- \Lx Q 3(x / )(G(x)) dx >\J ni Q2{x')(G(x) tan ) dx. 

Above, we used the fact that XhG h converges weakly in L 2 (Q 1 , M 3x3 ) to G (compare with the 
convergence in (5.2)) and the nonnegative definiteness of the quadratic forms Q 3 (x'), following 
from A{x') being the minimizer of the mapping W, as above. 



By (j5_3|: 

Q 2 {x')(g{x',x 3 ) 



I an 



Q 2 {x') (g{x\ 0) tan ) + 2x 3 C 2 {x') (G(x', 0) tan , [AQ^Vn] 
+ xlQ 2 {x')([AQ~ l Vn) tan {x') 



tan\X 



The second term above, expressed in terms of the bilinear operator C 2 {x') representing the qua- 
dratic form Q 2 (x'), integrates to on the domain fi 1 symmetric in x 3 . After dropping the first 
nonnegative term, we arrive at: 

liminf ^l\u h ) > ^ j^Q2{x')([AQ- l Vn] tan {x')) dx' . 
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This yields the formula in (iii), as AQ" 1 £ 50(3) so AQ- 1 = (QA- 1 ) 7 = A~ l Q T . The proof of 
Theorem |2.4|is now complete. ■ 



6. TWO LEMMAS ON THE QUADRATIC FORMS Q3 AND Q 2 

We now gather some facts regarding the quadratic forms Q3 and Q 2 . First, it is easy to 
notice that the tangent space to T{x) at A{x) coincides with the products of the skew-symmetric 
matrices with A(x), denoted by skew • A{x). Here 'skew' stands for the space of all skew matrices 
of appropriate dimension; in the present case 3x3. The orthogonal complement of this space 
equals to: 

= sym • A(x)~ 1 , 



T A {x)F{ x ) 

where 'sym' denotes the space of all symmetric matrices (again, of appropriate dimension). The 
quadratic, nonnegative definite form Q 3 {x') is strictly positive definite on the space above, and it 
depends only on the projection of its argument on this space: 

(6.1) Q3(x')(F) = Q 3 (x') {F {sym . A(xf) -i } F) . 

Lemma 6.1. We have: 

^{sym-A- 1 } 

with 63 = — /33 and: [ b\ b 
Proof. Since the projection 



F a p 
hi f'32 

Id + A 



/l3 
/23 
/33 



{sym.A;l} F "/3 

[bi b 2 ]A~l 

/23 ] — [ /31 /32 
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h 
bo 

b 3 
]A 



1 

a/3" 



is a linear operator, we will separately prove the above formula 



in two cases: when F a p = and when fij = 0. Notice first that 



{sym- A 1 } 



F = BA- 1 , for a 



symmetric matrix B, uniquely determined through the formula: 

VSesym = {F - BA' 1 ) : {SA- 1 ). 
Since the right hand side above equals to (FA^ 1 — BA~ 2 ) : S, we obtain: 



(6.2) FA 
Also, we notice the form of the matrix: 

(6.3) B 



BA 2 E skew. 



F>a/3 

bi b 2 



b 2 
h 



BA 



-1 



B <*P A af3 

[h b 2 ]A-J 
]A n8 . ThenF„ 



bi 
b 2 
h 



In the first case when f {j = 0, let B aj3 = [P {sym . A -i } F al3 ]A al 3. Then F a/3 A a p - B a/3 A a ^ £ skew, 
and the same matrix provides the only non-zero, principal 2x2 minor of the 3x3 matrix 
FA -1 — BA~ 2 , where B is taken so that all 6^ = and Bt an = B a p. By uniqueness of the 
symmetric matrix B satisfying (6.2), this proves the claim. 

In the second case when F a $ 



FA- 



BA- 



0, define B as in (6.3) with B a p 




0. The result follows, since: 



[ /31 /32 ] A 



aji 



fl3 
/23 
/33 



+ 



I [61 



A 



a/3 



bi 
In 



and (6.2) is equivalent to the conditions on bi given in the statement of the lemma. We remark 



that since A a ^ = [g a f3\ 1 is strictly positive definite, then the same is true for the matrix Id + j4 
which implies its invertibility. 



a/3' 
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Now, the quadratic and nonnegative definite form Q,2(x') is likewise strictly positive definite on 
the space sym • and: 

(6-4) Q2{x'){F tan ) = Q 2 {x') (P^ym.A^)- 1 }^) ■ 

Lemma 6.2. There exists linear maps b,c : M 2x2 — ► M 3 related by: 



c{F t 



tan ) 



and such that: 

(6.5) Q 2 (x'){F tan ) 



Q 3 (x') 



Id + A-J 




A 



a/3 



h 
b 2 
b 3 



■b(F tan ) 



Qsix') 



Ftan 





Proof. By (|6.1|), Lemma 6.1 and the definition of Q2 it follows that: 

'{sym- AllY 



Q2{x){F tan ) = min Q 3 {x 



-li Ftan 



h 
b 2 

?3 



Cl 

C>2 
C3 



dence we obtain t 



(6.1 ) and Lemma 



6.1 



provided that C3 = —63 and [ c\ C2 



[ [bl b 2 ]A~ a l 63 
;eco 

[b 1 b 2 I (Td + ,C1). which 



re fi rst equality in the representation (|6.5|). The second one follows again from 

2^ 

a/3 



is exactly the condition defining the vector c in the statement of the lemma. 



7. The recovery sequence - a proof of Theorem 2.5 



6.2 



wc 



Following the reasoning in step 1 of the proof of Theorem 2.4 we first notice that the matrix 
field Q whose columns are given by: 

Q(x') = [d 1 y(x'),d 2 y(x'),n(x')\ ef(x'). 

Hence in particular: QA~ l G SO(3). With the help of the above definition and Lemma 
put: 

(7.1) d(x') = Q(x')A-\x') • C (^(V2/) T Vn(x')) G L 2 (n,R 3 ). 
Let d h G W 1 ' 00 ^,^ 3 ) be such that: 

(7.2) lim d h = d in L 2 (n) and lim h\\d h \\ w i,oo = 0. 

h— >0 



limh\\d h \\ w i 

h— >0 



Note that a sequence c?' 1 with properties (7.2) can always be derived by reparametrizing (slowing 
down) a seq uenc e of smooth approximations of the given vector field d G L 2 (S7). 

Recalling (2.5), we now approximate y and n respectively by sequences y h G Pv a2 '°°(Q, IR 3 ) and 
fi A G ^^^(r^M 3 ) such that: 

lim \\y h - y\\ W 2,2(Q) = 0, lim \\n h - n\\ w i,2 {n) = 0, 

(7.3) ^(ll/llv7 2 >-(o) + ll^llwL-cn)) < eo, 

lim ^ \x' G ft; y^x') ^ y(x')} U |x' G ft; n^x') / ra(x')} = 0, 
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for a sufficiently small, fixed number eo > 0, to be chosen later. The existence of such approxima- 
tion follows by partition of unity and a truncation argument, as a special case of the Lusin-type 
result for Sobolev functions in |2U| (see also Proposition 2 in [B]). 
Define: 



(7.4) 



u h (x',x 3 ) = y n (x')+x 3 n n (x') + ^d n (x'). 



hi i 



h(J\ 



X 



Note that each map: $7 3 x' i— > dist(Vu fc (x'), J-(x')) vanishes on O/j and is Lipschitz in $7, with 
Lipschitz constant of order 0(1/ h). Here, we let: 







h = |x' € ft; y h (x') = y(x') and n h (x') = n(x')} . 



For any point x' G fi\fi/,, we also have dist (a/, 0^) < C|fi\0/j|. The proof of the latter statement 
is standard, see for example [T8], Lemma 6.1 for a similar argument. As a consequence, by (7.3) 
we obtain l//i 2 dist 2 (x', $7^) — > and hence: 

(7.5) dist(Vu'V), F(x')) < 0(l/h)dist(x' , « h ) = o(l). 

The gradient of the deformation u h is given by: 



Vu h (x',x 3 ) = Q h (x')+x 3 A h 2 (x') + X ^D h (x') 



where: 



Q h (x') = Q(x') infi fc , A£(x') = d 1 n h (x'),d 2 n h (x'),d h (x') 



lim An = A 2 



D r < 



din, d 2 n, d 
d 1 d h ,d 2 d h ,0 



in L 2 (Q), 



Note that by (7.5 ) and the local C 2 regularity of W, the quantity W(x, Vu h (x)) remains bounded 
upon choosing h and eo in (7.3) small enough. The convergence in (i) Theorem 2.4 follows 
immediately. 

We now prove (2.6). Using Taylor's expansion of W in a neighborhood of Q(x'), we obtain: 

2 



\\ W (x,Q(x') + hx A r 2 l (x') + h 2 ^D h (x')) dx + ^ [ W(x, Vu h (x)) dx 
^ ifli \ - • 2 ) h 2 J a i\ai 

= f (^V 2 ^(x / ,-)|Q(,o(^4(^)^34(^)) + ^(^)) da; + ^|fi\fi h |. 

Here the reminder lZ h converges, by (7.2 ), to pointwise almost everywhere, as h — > 0. Ther efore, 
recalling the boundedness of W(x, Vu h (x)) we deduce by dominated convergence and (7.3) that 
the above integral converges, as h — > 0, to: 

l - J x 2 V 2 W(x', ■) lQ(x , ) (A 2 (x'),A 2 (x')) ^=\J ni 4Qs(x') (AQ^Ao) dx 

= ki J n ^ x '){ A ~ l Q TA ^) dx ' = ^J n 22(^(Vy) T Vn) dx'. 
where we applied frame invariance, (7.1) and (6.5). ■ 
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8. Conditions for existence of W 2 ' 2 isometric immersions of Riemannian metrics - 

a proof of Theorem 12.61 



The assertions in (i) follow directly from Theorem 2.4 and Theorem 2.5 It remains to prove 
(ii), which clearly implies (hi). 
Assume that lim/j^o j^I h (u h 



by Theorem 



2.4 



for some sequence of deformations u h G W 1,2 (Q h , M 3 ). Then, 
there exists a metric realization y G W 2 ' 2 (fl, M 3 ) such that: 



where II = (Vy) T Vn is the second fundamental form of the image surface y{^l). Recalling (6.4) 
we obtain: 



= Q 2 (x / )(^n) =Q 2 (x' 



sym-A-^ Q/ 3 J 



i(Ajn)) Vx'eo. 



Since the quadratic form Q,2(x') is nondegenerate on sym • A a jtj, it follows that: 



a{3 



X.A-l(Kln)=Q, 



.1), so: 



for the symmetric matrix B G M 2x2 satisfying: 

(A^H - BA'l) : (5A-1) =0 VS 1 G sym. 

The above condition is equivalent to A~^I1A~^ — BA~^ G skew, but B = in view of 

Kl^Kl 6 skew. 

Since II G sym, there must be II = and therefore indeed effectively y : O — ► M 2 . 

On the other hand, if y G W 2 ' 2 (tt, M 2 ) is a 2d realization of [g Q( g] then clearly II = (Vy) T Vn = 0, 
so for the recovery sequence corresponding to y and constructed in Theorem |2.5[ there holds 
linnwo j?I h {u h ) = I(y) =0. ■ 
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